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Abstract 

Using the published Kamiokande binned data for sub-GeV and multi- 
GeV atmospheric neutrinos, we have searched for the optimum set of three 
flavor neutrino oscillation parameters within the constraints of reactor ex¬ 
periments. It is found that y 2 is minimized for (Am^, Am| 2 , #12, #13, 6*23) = 
(7.4xl0 _1 eV 2 , 2.6xl0“ 2 eV 2 , 2°, 3°, 45°) and (2.6xl(T 2 eV 2 , 7.4X10" 1 eV 2 , 
0°, 87°, 46°) with xLm — 87-8 (76%CL). The sets of parameters which are 
suggested by the two flavor analysis turn out to be close to the minimum 
(0.2a). 
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There has been much interest in atmospheric neutrinos 1,0. g g g. 
which might give evidence for neutrino oscillations. While NUSEX H and 
Frejus jSJ have reported consistency between the data and the predictions of 
the atmospheric neutrino flux |], ||, Kamiokande, IMB and Soudan-2 have 
reported a discrepancy. In particular, the Kamiokande group claimed that 
their data suggests that the mass squared difference is of order lCT 2 eV 2 [§]. 
The reason that they have obtained a narrow region for the mass squared 
difference is because they have used both the binned data in the sub-GeV 
and the multi-GeV energy regions. Namely, while the multi-GeV data show 
remarkable zenith angle dependence, the sub-GeV data have little zenith 
angle dependence. This indicates that the mass squared difference relevant 
to the neutrino oscillations in the atmospheric neutrinos cannot be much 
smaller than or much larger than 10 _2 eV 2 . The momentum spectra of the 
sub-GeV atmospheric neutrinos (Fig.l in [l|) also supports this argument. 

People have studied neutrino oscillations among three flavors 0- and 
it has been shown recently |TT that one can easily get strong constraints 
for the mass squared differences and the mixing angles if one assumes a 
mass hierarchy. The original analysis of the atmospheric neutrino data by 
the Kamiokande group II was based on the framework of neutrino os¬ 
cillation between two flavors and it is important to see what happens if we 
analyze the data in the three flavor framework. Much work has been done 
on the analysis of atmospheric neutrinos 0. 0. [14], [l5j| and among those 


that discuss quantitatively the binned data of the multi-GeV neutrinos are 
[L3|, [L4], [K|. However, these works did not take the sub-GeV binned data of 
Kamiokande [l| into account, and the allowed regions for the relevant mass 
squared difference obtained in there are wider than those in |J. In this pa- 
perQ we will perform a three flavor analysis of the published data |], [2] of the 


2 This paper supersedes jlii 
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Kamiokande atmospheric neutrino experiment, along with the data of the 
reactor neutrino experiments [16|, [17]] . We will make full use of both the sub- 
GeV and the multi-GeV binned data as well as the data of reactor neutrino 
experiments [| 


171, and as we will see below, we have strong constraints on 


the mass squared differences and some of the mixing angles. We will take 
the smearing effects into consideration, and evaluate the number of events by 
summing over the energy and the zenith angle of neutrinos, to reproduce the 
original analysis by the Kamiokande group as much as possible. Throughout 
this paper we will restrict our discussions only to the data by the Kamiokande 
group |], 0, not only because the Monte Carlo result for the neutrino energy 
spectrum is available only in Ref. i0. but also because this is the only 
data which gives both the upper and the lower bound on the mass squared 
difference of neutrinos. 

We start with the Dirac equation for three flavors of neutrinos with mass 


in matter [|18 
. d 

i— \k(x) = 
dx 


Kdiag (o, Am^/217, Aml 1 /2E^j U 1 + diag (H(x), 0, 0) T(x), (1) 


where we have taken the ultra-relativistic limit and have subtracted a term 
proportional to the unit matrix, Araf - = m 2 — m 2 is the mass squared differ¬ 
ence of the neutrinos with energy E , T(x) = (z/ e (x), ^(x), u t (x)) t is the wave 
function of the neutrinos in the flavor basis, and A(x) = \Z2GpN e (x) stands 
for the effect due to the charged current interactions between v e and electrons 


in matter [|18[]. Here, 

/ 


U = (U aj ) = 


C 12 D 3 

S12C23 — C12S23S13 
V S12S23 — C12C23S13 


S12C13 S13 

D2C23 — S12S23S13 S23D3 

-C12-S23 — S12C23S13 C23C13 


( 2 ) 


with Cij = cos Oij, Sij = sin % (a = e, p, r; i,j = 1,2,3) is the orthogonal 
mixing matrix of neutrinos. For simplicity, we will not discuss the CP vio- 
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lating phase of the mixing matrix here. Q. Without loss of generality we can 
assume m\ < < m§, so that we have > max(Am 21 , Am| 2 ), and we 

will use Am 21 , Am| 2 > 0 as two independent parameters. 

The number of the expected charged leptons i a ( £ Q = e or p ) with 
energy q from a scattering u a N —> £ a N' ( a = e, /i) is given by 

r°o rQmax r 1 /• 1 r 2 ir 

N(£ a ) — tit V / dE / dq d cos© / dcos 6 * / dy? e Q (g) 
jo Jo J-i J-i jo 

P— 

d 3 F /3 (E, 9) _ d 2 a a (E, q) d cos^A ■ F 

dE d cos 9 dip dq dcos'ijj d cos 9 ^ ’ 

Here d 3 Fp/dEd cos 6dp is the flux of atmospheric neutrinos z/g with energy E 
from the zenith angle d , tit is the effective number of target nucleons, e Q (g) 
is the detection efficiency function for charged leptons £ a , d,a a /dqd cos is 
the differential cross section of the interaction z/ a IV —> AW' (a = e or /i), 
and 0 is the zenith angle of the direction from which the charged lepton £ a 
comes (See Fig. 1). 

(Insert Fig.l here.) 


P{yp —> zz Q ; E , 6) is the probability of z/^ —> u a transitions with energy E 
after traveling a distance L = \J(R + h) 2 — R 2 sin 2 9 — R cos 9, where R is 
the radius of the Earth, h ~15Km is the altitude at which atmospheric 
neutrinos are produced. 

As for the analysis of the sub-GeV data, almost all the information which 
is necessary to get the right hand side of (||) is available in the published 
references. We have used the differential cross section da a /dqd cos i/j in [|H| , 
and the detection efficiency function e a (q) is given in |2(J. The ffux of atmo¬ 
spheric neutrinos d 3 Fp/dEd cos 9dp without the geomagnetic effects is given 


3 Even if we include the CP violating phase 6 of the mixing matrix, the effect of 6 
always appears in the combination of C13 sin 2013 shA. sin2$i3 has to be small because of 
the constraints from the reactor experiments. 
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in [||] but we have used the flux which has been obtained with the geomag¬ 
netic effects [Q. There are two important types of binned data in [lj from 
which we derive the dependence of neutrino oscillations on the energy and the 
path length of neutrinos: Fig.l (a) and (b) in [[1] give the binned data with 
respect to energy of the outgoing charged leptons, while Fig.3 in |T| gives the 
binned data (the ratio (/i/e)d a ta/(/V e )Mc) with respect to the zenith angle. 
These binned data are expressed as 


w , rcos@j +1 roc rq m ax r 1 /-2 n 

Y a = (1 + a) ht V) Cn / d cos© / dE / dq dcosO / dp 

J , ./cose, Jo Jo J -1 Jo 

. d 3 Fs(E,9 ) d 2 a a (E, q) dcos-0 ... 

(q) ' jjp \ - , 1 -pr P(yp -► E, 9) 4 

dE d cost/ dp dq dcosip d cost) 


= (1 + a) n T £ Cp I ‘ dq I dE 

ka + l JO J-l 


P=e,V 


d cos© / dcos9 / dp 
J -i Jo 


xe a (q) 


d 3 F/ 3 (E,9 ) d 2 a a (E,q ) dcos0 


dE dcos9 dp dq dcos'ijj d cos© 


P(up^u a] E,9), (5) 


respectively, where (1 + a) is a factor for the uncertainty in the absolute 
normalization, C^ = 1+(3/2, C e = 1 — /?/2 for the uncertainty in the relative 
normalization of the and v e flux, and cos Qj = (2 j — 7)/5 (1 < j < 5), 
q a = (a + 1)/10 GeV (1 < a < 10). 

To reproduce the analysis of the multi-GeV data by the Kamiokande 
group, one needs the detection efficiency function e a (q), which is not given 
in |Q]. However, the quantity 

Mmsx f 1 r 1 r2ir 

g a {E) = tit / dq e(q) / d cos© / dcos9 / dyj 
jo J-i J-i Jo 

^ d 3 F a (E,9 ) d 2 a a (E,q) dcos'ijj 
dE d cos 9 dq d cos 0 d cos 0 


is given in Fig.2 (d)-(f) in Ref. The zenith angle dependence np(E,Q) 
of the atmospheric neutrino flux for various neutrino energy E has been 
given in Ref. || in detail. Here we multiply the quantity g a (E ) by the 
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zenith angle dependence n a (E,Q) in Ref. || || with suitable normalization 
N' = 1/ fgdcosO n a (E,0), and insert the Gaussian factor with la = 17° 
to simulate a smearing effect of the detecting neutrinos. Thus we adopt the 
following quantity for the multi-GeV analysis. 

X? = NN\l + a) Y, C pj~ dcosQ J dE J J dcos9dtp 

(—tan 2 -0/2/tan 2 ^0/2) ( ^ 

! 9 - P{vp -»• E, 9), (7) 

(1 + cos -0) \J sin 0 — sin 9 sin p> 


0=e,fi 

9 p(E)np(E, 0) exp 


where D is the region of ( 9 , ip) in which the argument of y sin 2 0 — sin 2 9 sin 2 <p 
becomes positive, N = 2/ v / 7rtan('0 o /2) is the normalization factor such that 


[ * 1 exp(— tan 2 ip/2/ tan 2 ^o/2) 

iv / dcosw- --7777-7-7-7- = 

J- 1 (1 + COS'?/’) 3 / 2 (1 — COS '^) 1 / 2 

and we put ijj 0 = 17°. The relation between ip and 0 


( 8 ) 


cos ip = 


cos 9 cos 0 — sin 9 cos p>\j sin 2 0 — sin 2 9 sin 2 ip 


(9) 


1 — sin 2 9 sin 2 <p 

can be obtained by spherical trigonometry, and the measure d cos pj/d cos 0 
gives factors in the numerator in (|7|). (0) is not exactly the same quantity as 
the one in the original analysis ||, but this is almost the best which can be 
done with the published data in 0 . 

Several groups [[7], |j have given predictions on the flux of atmospheric neu¬ 
trinos but they differ from one another in the magnitudes, and the Kamiokande 
group assumed that the errors of the overall normalization 1 + a and the rel¬ 
ative normalization 1+/3/2 are cTq,= 30% and 073=12%, respectively. Thus we 
define the total y 2 as plj 


a 


P 


X o ^ 2 Xsub—GeV "b Xmulti—GeV f XBugey ~b XKrasnoyarsk; 


cr. 


a 


( 10 ) 
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where 


10 

E 

a=e,ii a = 1 


yet 

1 a 


Xsub-GeV = 2 E E [ Y a * Va ~ Va ln + 


Va 


5 1 
3 =1 3 


'Xi 7 # ; 

,e 1 


^ 7 e /2/f 


Xmulti—GeV 


Xb 


ugey 


( 11 ) 
( 12 ) 

a=e,nj =1 V +7 / 

is given by (9) in |T(J, and ^Krasnoyarsk i s defined in terms of the eight 


= 2 E e(e 


data and the three parameters N\, N 3y Nb defined in [[L7 


Here x“, yf 


a = 


e, y;l < j < 5) are the mnlti-GeV and the sub-GeV data for each zenith angle 
cos0j < cos© < cos0j+i, respectively, and (a = e,p; 1 < a < 10) are 
the sub-GeV data of the charged leptons with the energy q a < q < q a+ \. It is 
understood that a 2 /<r 2 + /? 2 /cx| + Xsub-GeV+Xmuiti-GeV is ^ ie optimum value 
with respect to a, /3, y| with respect to the parameters A, 6, aj(j = 1, 2, 3) 


defined in H, and XKrasnoyarsk with respect to the parameters Ni,N 3 ,Nb. 
The theoretical predictions V“, Yj* (1 < j < 5), V a a (1 < a < 10), (a = e, //) 
depend on five free parameters (Am| 1; Am| 2 ; $ 12 , # 13 , $ 23 ), where Arn l3 = 
rrif— m 2 , so ( p~Q| ) is expected to obey a y 2 distribution with 103—5=98 degrees 
of freedom (10x2+5=25 for sub-GeV, 5x2=10 for multi-GeV, 60 for Bugey, 
and 8 for Krasnoyarsk). The number of degrees of freedom in the present 
analysis for the sub-GeV and the multi-GeV data is smaller than in the 
original one by the Kamiokande group (10x11=110 for sub-GeV, 5x8+5=85 
for multi-GeV). 

The value of y 2 is affected to some extent by the presence of matter, and 
it is necessary to take into consideration the contribution of the second term 
in ([!]). In the case of the multi-GeV analysis we have solved (jl]) numerically 
for each E (0.9 GeV < E < 100 GeV) and evaluated the number of events for 
a given range of the zenith angle. As for the sub-GeV case, however, it turned 
out that the second term in (jl) is negligible for the entire energy region of 
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the sub-GeV neutrinos with Am?- under consideration (A(x) -C Am? /E), 
so we have ignored the matter contribution in the sub-GeV analysis. We 
have also taken into account the effects of the particle misidentihcation error 
and the backgrounds (Y e —» Y e + 0.04Y M , Y 11 —■> 0.96W for sub-GeV [||], 
X e —■> X e + 0.08A"G X' Ji —■> 0.92A"^ for multi-GeV |22|] 0 ). We have almost 
reproduced the Monte Carlo results by the Kamiokande group, such as the 
energy spectrum of the sub-GeV (Fig. 1 in |l]]) the zenith angle distributions 
of the sub-GeV (Fig. 2 in [(TJ) and the multi-GeV data (Fig. 3 in |2[]). 

We have meshed each Am? region into 18 points (Am? = io( 2 ^- 24 )/ 5 , (0 < 
£ < 17)) and each Q %] region into 4 points % = £ 77 / 6 , (0 < £ < 3), and 
the evaluated the value of y 2 . Furthermore, using the grid-search and the 
gradient-search methods described in Ref. [^4|], we have found that y 2 has 
the minimum value for 


(Am^, Am^) = (7.4 x 10 1 eV 2 ,2.6 x 10 2 eV 2 ) 

(^ 12 ,^ 135 ^ 23 ) — (2°, 3°, 45°) 

(a, (3) = (3.1 x 10~\ -6.5 x 10~ 2 ) (13) 

with y 2 lin = 87.8 and 

(Am 2 !, Am 2 2 ) = (2.6 x 10~ 2 eV 2 ,7.1 x 10~ 2 eV 2 ) 

(012,013,023) = (0°, 87°, 46°) 

(a, p) = (3.1 x 10 _1 , -6.9 x 10~ 2 ) (14) 

With Xmin = 87 - 8 - 

The number of degrees of freedom of our analysis is 98, so the value of the 

reduced chi square is 0.9, which corresponds to 76 % confidence level. This 

suggests that our fit in the present analysis is good. The 68 % CL, 90 % CL 

4 The 4% uncertainty (out of 8%) of the multi-GeV data is due to the /x-like neutral 

current events pi|. 
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allowed regions of the parameters Am| 1; Am § 2 with unconstrained $ 12 , # 13 , 6*23 
are given by y 2 < y 2 lin + 5.9, y 2 < y 2 nin + 9.2 respectively, and are given 
in Fig. 2, where the asterisks stand for the sets of the parameters for xL\ n - 
The reason that we have obtained a rather narrow region for (Am|i, Am§ 2 ) 
is because we have taken into consideration the multi-GeV binned data, the 
sub-GeV binned data and the reactor data together: The significant zenith 
angle dependence of the multi-GeV data gives certain lower and upper bounds 
on Am-j and the little zenith angle dependence of the sub-GeV data gives 
more constraint on the lower bound. In general X^uiti-GeV prefers a large 
value of sin 2 26*13 (if Am 21 < Am 2 2 ) or a large value of 4|[/ e i| J (l — \U e \\ 2 ) 
(if A 77 ^ 2 ! > A'/ 7 i 2 2 ), but the reactor experiments exclude a region with large 
value of either factor if the larger Ara| is of order lCG 2 eV 2 . 

(Insert Fig.2 here.) 

In Fig. 3 we have plotted y 2 as a function of 0 13 , putting Am 21 = 0 with 
6*23 unconstrained. This is a situation which is realized by the two flavor 
analysis of the solar neutrino problem [|l^, Irrespective of whether we 
consider the vacuum solution (Am 21 ~ (9(10 _ 11 eV 2 )) or the MSW solution 
(Am 21 r-u 0( 10 5 eV 2 )) for the solar neutrino, the mass squared difference 
Ai s negligible compared to the contribution of the matter effect A(x ) in 
(|ip and the other mass squared difference A m 2 2 , which should be at least of 
order 10^ 2 eV 2 to account for the zenith angle dependence of the multi-GeV 
data. In this case 0 12 does not appear in y 2 . In particular, for 0 13 = 0 we 
have x 2 ~ xL = 2 - 2 , so we conclude that this set of parameters falls within 
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0.2cr for all three types of the solar neutrino solutions [[18], ^5| []. 


(Insert Fig.3 here.) 


In this paper we have analyzed the sub-GeV and the multi-GeV atmo¬ 
spheric neutrino data by the Kamiokande group based on the framework 
of three flavor neutrino oscillations along with constraints from the reactor 
experiments, and have shown that at least one of Amf- should be of order 
10~ 2 eV 2 at 90 % confidence level. We have also shown that the popular 
set of parameters ( (Am^, sin 2 2 #i 2 ) = (Am 2 , sin 2 29) Q , (Am| 2 , sin 2 2023) = 
(Am 2 , sin 2 20) atm , d 13 = 0) fall within 0.2cr, which is very close to the best 
fit. The minimum value of y 2 is 87.8 for 98 degrees of freedom, and the 
fit based on the hypothesis of neutrino oscillations is good. If we combine 
the present results with the solar neutrino experiments, then we get even 


stronger constraints, which will be reported somewhere [26 
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Figures 

Fig.l The parametrization of angles in the interaction u a + N —> i a + N'. 

Fig.2 The allowed region for Am | l5 Am | 2 with unconstrained. The solid 

and dashed lines stand for the regions at 90 % and 68 % confidence 
level, respectively. The asterisks stand for the sets of the parameters 
for the best fit. 

Fig. 3 x 2 as a function of d 13 with Am^ = 0 and d 2 3 unconstrained, y 2 is 
independent of 9 12 in this limit. 
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